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Abstract 

In a recent paper, Li and Melrose have claimed that the splitting — due 

to relativistic effects — of the Alfven surface in an axisymmetric pulsar wind 

y—( does not occur. Here we refute this claim by showing that, unless the solution 

^ that describes the flow along each open magnetic fleld line passes through 

Oh 

< 



the pure Alfvenic point (which is one of the manifestations of the splitting 
of the Alfvenic point), it would not be physically viable both at the surface 
of the star and at inflnity. 
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O !• Introduction 

Under the conditions of axial symmetry and time-independence, several 
flrst integrals of the system of magneto-hydrodynamic equations that govern 
a relativistic pulsar wind are known (Ardavan 1976). The conservation laws 
2 for mass, energy and angular momentum, together with Ferraro's isorotation 

law, enable us to establish a relationship between the values of any two flow 
^ variables along a given open magnetic fleld line. Each of these relationships 

can generally be written in the form of an algebraic equation for one variable 
as a function of the other. For instance, the Lorentz factor of the plasma 
satisfles a quadratic equation whose coefflcients depend on the plasma density 
and the distance of the fleld point from the axis of symmetry. Amongst the 
various solutions of such an algebraic equation, there is normally none that 
can describe the entire flow on its own: the solution that satisfles the required 
boundary conditions at the surface of the star fails to have an acceptable 
asymptotic behaviour at inflnity, and the solution that behaves properly at 
inflnity becomes unphysical at the star. 
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A physically acceptable description of the flow can therefore be obtained 
only by matching two (or more) solutions of the algebraic equation in question 
along a given magnetic fleld line. If the flow is to be free of discontinuities, 
however, these distinct solutions must be matched across surfaces at which 
the discriminant of the equation vanishes without changing sign. Such zeros 
of the discriminant of each equation deflne the loci of the critical (or singular) 
points through which the composite solution describing the actual flow should 
pass. The Alfvenic singularity — at which the poloidal component of the flow 
velocity equals the local propagation speed of Alfven waves — is one of these 
points. 

Of the several critical points thus arising from the various algebraic 
equations that are implied by the flrst integrals of the equations of motion, 
there are two — the so-called pure Alfvenic and intermediate critical points — 
which coalesce and coincide with the Alfvenic point in the non-relativistic 
regime. Not only these offshoots of the Alfvenic point, but every one of the 
critical points that arise in the context of constructing physically acceptable 
composite solutions must necessarily occur within the flow. 

There are other contexts in which the Alfvenic point and the manifesta- 
tions of its relativistic splitting arise. From the flrst integrals of the equations 
of motion, one can also derive expressions for certain flow variables which 
have the form of a ratio. It so happens that the zeros of the denominators 
in such ratios often coincide with the critical points of the flow as identifled 
above. Li & Melrose (1994) deflne the critical points or singularities of the 
flow by means of such ratios alone and, from the fact that the numerators 
and the denominators in these ratios have common factors when expressed 
in terms of certain variables, draw the conclusion that the pure Alfvenic and 
the intermediate singularities are spurious. 

If the only argument for the occurrence of these critical points were 
the possibility of the vanishing of the denominators in the expressions that 
one obtains for various flow variables, then the observations made by Li and 
Melrose would have implied that this argument does not suflice to establish 
the occurrence of the pure Alfvenic and the intermediate critical points: the 
zeros of the numerators and the denominators in the expressions in question 
may automatically coincide. However, as we shall see in §3, the main argu- 
ment for the splitting of the Alfven surface is based on the global behaviour 
of the solutions of certain quadratic equations and is in fact totally indepen- 
dent of the ratios by means of which Li and Melrose deflne the singularities 
of the flow. 

2. Formulation of the problem: the Alfven surface 
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The magneto-hydrodynamic equations which govern a relativistic pul- 
sar wind have the foUowing first integrals when the flow is time-independent 
and axially symmetric: 

V = K.B + rtoe^, (1) 
= F, (2) 
—rB^/{4:Tr) + r^fpv^K = G, (3) 

.Jl_!2^!!!^)=expH (4) 
c c 

(see Ardavan 1976). In these expressions, v is the flow velocity, B is the 
magnetic fleld, p is the rest-mass density, 

l = {l-vyc')-h (5) 

is the Lorentz factor, k is a scalar function of position, r is the distance from 
the axis of symmetry, e^p is the unit vector associated with the azimuthal 
direction and the subscript (f designates the toroidal component of a vector. 
The quantities u^F^G, and H remain constant along the magnetic lines of 
force. If we exclude the possibility of the differential rotation of the magnetic 
fleld lines, then u represents the constant value of the angular velocity of 
rotation of the rigid surface of the central neutron star. The remaining 
constants G and H are related to the rates of transport of mass, angular 
momentum and energy per unit flux-tube, respectively. 

Equations (l)-(5) comprise six scalar relationships between the values 
of the coordinate r and the seven variables Vp^v^p, Bp, B^p, n and 7 along 
a given magnetic fleld line [vp and Bp denote the poloidal components of v 
and B). They can be employed to express six of these variables as functions 
of the seventh variable and r. In particular, one can express 7 as a function 
of AttF^k, = A'lT^pVp'^ / Bp and f = rto / c to obtain 



[l — f\ — AttF^k) exp H 
(1 - f^)(l - P - AttF-ik) 



7 = .r zr:: (6) 



in which 

r^ = (l + ^j^) • (7) 

The covariant representation of the statement that the velocity of the plasma 
in the polidal direction equals the Alfven speed in the present case reduces 
to 

4Tr-fpVpyBp^ = l-f^ (8) 
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(see the Appendix in Ardavan 1976). At the Alfven surface defined by (8) 
the denominator in (6) vanishes. Requiring that 7 should remain finite at 
this surface, i.e. that the numerator of the expression in (6) should vanish 
simultaneously with its denominator, we find that the open lines of force of 
the magnetic field cross the Alfven surface at r = f^. 

The Alfven cylinder (8) (whose generator is not in general a straight 
line) constitutes the surface of parabolic degeneracy of the system of par- 
tial differentail equations that governs the axisymmetric pulsar wind: this 
system, which undergoes a change in type from elliptic to hyperbolic across 
the magnetosonic surface Air^pVp'^ + Bp — = 0, becomes parabolic also 
at r = r^, within its domain of ellipticity (see Ardavan 1979). So, another 
way of inferring the occurrence of the Alfvenic singularity is to analyze the 
characteristics of the governing field equations. The calculation presented in 
§4 of Li & Melrose (1994) is part of the linearized version of such an analysis. 

The elimination of B^p between (1) and (3) results in an expression for 
v^p which has a structure similar to that of the expression on the right-hand 
side of (6): 

ru — AnGn/r 



(9) 



In the non-relativistic regime, where the denominators in (6) and (9) coincide, 
this expression, too, becomes singular at the Alfven surface and has to be 
regularized. 

In the relativistic regime, however, the surface 

1 -AirF-fK = (10) 

at which the denominator in (9) vanishes is distinct from the Alfven surface 
and the question arises as to whether this surface ever occurs within the fiow. 
If the so-called pure Alfvenic surface that is defined by (10) exists, then at 
this surface we should also have 

ru -AnGK/r = 0, (11) 

for the toroidal component of the fiow velocity, should be everywhere 
finite. 



3. Occurrence of the pure Alfvenic singularity 

That the fiow does in fact pass through surface (10) can be seen by 
considering the relationship between the Lorentz factor 7 and the density p 
of the plasma along an open field line. We show in this section that 7 is given 
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as a function of p and r by a quadratic equation neither of whose solutions 
provides a physicaUy acceptable description of the entire flow on its own. 
The solution that describes the flow close to the star is different from the 
one which has the correct asymptotic behaviour at inflnity so that the two 
solutions of this quadratic equation have to be matched across the surface 
at which the discriminant of the equation vanishes. The surface at which 
the discriminant of the quadratic equation in question vanishes, on the other 
hand, is precisely the same as the pure Alfvenic surface that we encountered 
in (10). 

If we use (1) and (4) to eliminate B^p and v^p from (3), and express n 
in the resulting expression in terms of /?, we arrive at the following quadratic 
equation for 7: 

-[{l-f^)R + {l-f\)-^]r + R = Q, (12) 

where 

r^e-% (13) 
i? = (47rF2e^)-V, (14) 
and rA is deflned in (7). Of the two solutions 

r± = i[(l - P)R + (1 - f\)-^] ± {i[(l - P)R + (1 - f\)-^f - rY^ (15) 

of this quadratic equation, only r_ passes through the Alfvenic singularity: 
when r = r^, it is r_ that equals (1 — f\)R and so satisfles the Alfvenic 
condition (8). 

The solution r_ , however, does not yield a positive value for the Lorentz 
factor at inflnity: it either vanishes or becomes negative as r — > 00. Since 
the plasma density p — and hence R — reduces to zero as r — > 00, (15) yields 
lim^^oo r_ = if the limiting value of —f^R + (1 — f\)~^ is non-negative 
and yields lim^^oo T- < if the limiting value of —f^R + (1 — f\)~^ is 
negative. Not only does it not yield a physically acceptable value for the 
Lorentz factor at inflnity, but the solution r_ also predicts a flux of plasma 
angular momentum in r > 1 that is directed towards the star. 

To see this, let us note that v^p is positive, and so the flux of plasma 
angular momentum is outward (cf. Ardavan 1976), only if F > 1 [see (4)]. 
In order that r_ exceeds unity, i.e. that 

(1 - f^)i? - (1 - 2ri)(i -fi)-'>{[{i-f')R+{i- fir^r -mK m 

we need to have both 

(^l_f')R-{l-2f\){l-f\)-' >0, (17) 
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and 

f^R-f\{l-f^^)-^ >0, (18) 

where the second inequality follows from squaring (16). Inequality (17) can 
be satisfied in f > 1 only if | < < 1. In this regime, however, it is not 
possible to satisfy (17) and (18) simultaneously: the upper limit set on R by 
(17) is smaller than the lower limit on this quantity that is required by (18). 

On the other hand, the solution r_|_ which has the correct behaviour at 
infinity does not yield an outward flux of electromagnetic angular momentum 
in r < rA- For this flux to be directed away from the star, B^p should be 
negative (see Ardavan 1976) and so, according to (1) and (4), (1 — f^)r should 
be smaller than unity. In order that (1 — r^)r_|_ is smaller than unity, i.e. 
that 

{[{l-r^)R+{l-r\)-^f-mY^ < 2{l - )-^- {I -r\)-^- {I -r^)R, (19) 
we need to have both 

2(1 - r^)-! _ (1 _ f\)-^ - (1 - f^)R > 0, (20) 

and 

f > f-A, (21) 

where (21) follows from squaring (19). Inequality (21) cannot be satisfled in 
r < so that the solution r_|_, too, is incapable of describing the entire flow 
on its own. This solution fails both to pass through the Alfvenic singularity 
and to describe the braking effect of the magnetic fleld in the sub- Alfvenic 
region. 

A physically viable description of the flow thus entails the matching 
of r_ and r_|_ across a surface in r > r^. The two solutions r_ and r_|_ 
of the quadratic equation (12) are equal to one another at points where the 
discriminant of this equation is zero. For the resulting composite solution to 
be free from discontinuities, it is therefore necessary that the two solutions 
match across a surface at which the discriminant of (12) vanishes, i.e. at 
which 

[(l-r')i? + (l-ri)-^]'-4i? = 0, (22) 

and (12) reduces to 

V = \[{l-P)R+{l-f\)-^]. (23) 

Equations (22) and (23) jointly yield the relationship V = R which deflnes 
the pure Alfvenic surface [see (10), (13) and (14)]. 
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The position of the pure Alfvenic surface now foUows from (11), (22) 
and (23): 



Equation (24) shows that the pure Alfvenic singularity always occurs beyond 
the Alfvenic singularity and that at this surface T = R = 1 [see (22) and 
(23)]. It also implies that, in the ultra-relativistic or massless limit where 
— > 1 [see (7)], the pure Alfvenic singularity is relegated to infinity (cf. 
Okamoto 1978). 

Note that the discriminant of equation (12) should vanish at the pure 
Alfvenic surface without changing sign, for otherwise (12) will not have real 
solutions on both sides of this surface. In other words, the discriminant in 
question and its gradient along a magnetic field line should vanish simulta- 
neously. This requirement further constrains the flow by specifying the value 
of B • V-R at the Alfvenic surface. It also ensures that r_ and r_|_ have flnite 
gradients at the pure Alfvenic surface: the numerators and the denominators 
in the expressions for B • VF-i- vanish simultaneously once this requirement 
is met [see (15)]. 

The above argument shows that, unless the solution that describes the 
flow along each open magnetic fleld line passes through the pure Alfvenic 
point, it would not be physically viable both at the surface of the star and at 
inflnity. This argument does not, however, guarantee that the solution that 
passes through the pure Alfvenic point is acceptable. A physically acceptable 
solution needs to meet not only the constraints associated with the pure 
Alfvenic point but also those associated with the other critical points of the 
flow. It turns out that the conservation laws (l)-(4) and the equation of 
the mixed type that governs the poloidal magnetic fleld do not admit of a 
continuous solution which can satisfy all such constraints (Ardavan 1979, 
Bogovalov 1994). 




(24) 
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4. Concluding remarks 

The question on which the conclusions of Li & Melrose (1994) disagree 
with those of Okamoto (1978) and Ardavan (1979) is whether the flow passes 
through the pure Alfvenic singularity deflned by (10) or not, i.e. whether 
there exists a surface within the pulsar wind at which the flow variables 
satisfy the relationships (10) and (11) or not. The answer to this question 
is independent of how one deflnes a singularity. Either the flow is further 
accelerated after passing through the Alfvenic point (8) to attain the higher 
value of the Lorentz factor that is required by the pure Alfvenic condition 
(10) or it is not. The answer to the question is dictated only by the physical 
conditions at the surface of the star and at inflnity and by the requirements 
of the conservation laws (l)-(4). 

What Li and Melrose observe is that the form assumed by expression 
(9) does not by itself establish the occurrence of the pure Alfvenic singularity 
because the numerator and the denominator in this expression may have a 
common factor. Even if the only argument for the occurrence of the pure 
Alfvenic point had been the possibility of the vanishing of the denominator 
in (9), this observation would not have ruled out the occurrence of the pure 
Alfvenic point as Li and Melrose claim; it would have only rendered the 
argument in question inconclusive. As it is, the possibility of the vanishing 
of the denominator in expression (9) is not in fact the only argument for the 
occurrence of the pure Alfvenic point. 

One can establish the occurrence of the pure Alfvenic point indepen- 
dently of expression (9) and only use this expression to further constrain 
the flow. The constraint that follows from requiring the numerator in (9) to 
vanish simultaneously with its denominator, moreover, has as much physical 
signiflcance as that which has to be satisfled at the Alfvenic point: T = R = 1 
and (24) relate the constants of motion (i^, G, H and uj) to the values of the 
flow variables at the singular point in exactly the same way as do (8) and 
f = va- 

The flaw in the reasoning of Li & Melrose does not lie solely in their 
assuming that the splitting of the Alfven surface has to be inferred from the 
ratios they consider. They also overlook the fact that different singulari- 
ties of the flow are exhibited by the relationships between different pairs of 
flow variables (Ardavan 1979). That the numerator and the denominator in 
the expression for v^p have common factors when expressed as functions of 
the coordinate x (= rto / c) and the flow variable y (= AirF'^^/p) does not 
imply that the pure Alfvenic singularity is spurious; it only means that the 
relationship between the particular flow variables v^p and y is singularity-free. 

The isolated Alfvenic singularity of the non-relativistic stellar wind in 
the present case splits into three distinct singular points: the pulsar wind also 
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includes an intermediate singular point, lying between the Alfvenic and the 
pure Alfvenic singular points, at which Air^pK^ = 1 — fv^/c. The analysis 
establishing the occurrence of this intermediate singular point entails the 
study of the global behaviour of the solutions of a quadratic equation which 
relates to p. This analysis closely parallels that presented here and has 
already been outlined in Ardavan (1979). 
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